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ABSTRACT 


This thesis investigates the synthesis of smoothed impulse responses using sampled 
data of the input and output of random noise driven electromagnetic systems. Special 
interactive software was developed for NPS’s time domain electromagnetic scattering 
laboratorv. The system performs signal acquisition, synthesizes time and frequency do- 
main scattering responses using broad band random noise and provides results as easily 
evaluated graphic displavs. Attempted validations of the svstem are made bv comparing 
synthesized impulse responses for mucrowave filters and transient scatterers to alternate 


experimental and computational data. 
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I. INTRODUCTION 


A. OVERVIEW 

The objective of this research is to demonstrate the viability of performing high- 
resolution impulse response scattering measurements using a broad-band noise source. 
The development of laboratorv facilities for high-resolution impulse response scattering 
measurements has generally proceeded using two major techniques. The first of these 
emplovs a stepped-frequency coherent oscillator and a vector (magnitude and phase) 
receiver. Impulse response target characteristics are obtained via inverse Fourier trans- 
formation of the frequency domain data. The second method obtains inipulse response 
measurements directly in the tune-domain by use of a repetitive fast-pulse target jllumu- 
nation with a sampling-scope acting as the receiver. The time-domain approach for 
Scultering imeasurements offers a viable alternative to the more prevalent conunuous 
Wave approach. Transient scattering measurements provide wavelorms that can be more 
M@imectiv interpreted as to cause and eflect and allow cxact range-gating of target re- 
sponses for climmation of unwanted clutter (Ref. 1. 2]. 

The development of digital samnpling oscilloscopes and broad band noise generators 
made possible a third method for Ingh-resolution impulse response scattering measure- 
ments, This technique ts based on the estimation of the crosscorrclution function be- 
tween input and output using a broad-band noise source as the transmitter. Thus 
methodology has been employed in measuring the mechanical impulse response of large 
structures, such as bridges and burdings [Ref. 3]. 

The practical advantages of this third technique are two-lold. The first is the wide 
baudwidth and high power available from noise sources which are available at much 
lower cost than comparable stepped-frequency sources. A conclusive demonstration of 
the viability of noise source impulse response measurements mav lead to further refine- 
ments and ultimately to commercial marketing of this technology. The second advan- 
tage 1s related to the use of notse-source illumination for tactical and strategic radar 
applications. An obvious benefit would be the masking of the radar interrogation signal; 
this would appear at the target as either jamming or interference. Confusion would re- 
sult as to proper countermeasures to be emploved by the target. Furthermore. stnce this 
method uses random noise, which is orthogonal to other signals and other icise 


(producing Zero crosscorrelation with these), a high SNR requirement could be achieved. 


This thesis is a continuation in a series of efforts in transient electromagnetic scat- 
tering that began in 1979 at the Naval Postgraduate School (NPS). The development 
of digital sampling techniques allowed the developinent of transient scattering ranges in 
the late 1960's [Ref. 4]. A transient scattering range having sullicient bandwidth and 
signal to noise ratio (SNR) to support radar target identification was untially con- 
structed at NPS in 1980 using a ground plane configuration [Ref. 5,6]. A very wideband 
(short pulse) free-field scattering range was then constructed in 1985 for implementation 
of target identification based upon natural resonances, as introduced by Mittra and Van 
Blaricum (Ref. 7]. This free field transrent range was validated by Manateciiae 


Wee ainiels| ere, 0). 


B. NOISE SOURCE IMPULSE RESPONSE NIEASURERIENT CONG@iREES 

In the noise source unpulse response measurement technique, a broad band noise 
signal is used as the input. The system impulse response is derived from measuring the 
crosscorrelation function between the input and output. 

The configuration shown tn Figure la has been used to unpiement the poise source 
nnpulse response measurement techmique m the NPS scattering laboratory. It was 
thought at first that the Tektronix 7804 Dreital Processing Oscilloscopes (DPQ) would 
allow stnultaneous measurements of both channels. After much effort. this was found 
to be untrue, as will be shown. AV new method was developed, using the measurement 
of the difference of the input and output time sequence. A digtial estimator of the 
crosscorrelation function 1s formed in the data processor for the measurement and the 
expectation of the crosscorrelation function, &,,(7), 1s obtained by further averaging of 
this estimator using N-block averages of the stored data. 

The complete scattering measurement using a noaise-source requires three steps. 
Looking at the scattering range svstem representation in Figure [b, the target transter 
function //(f) will be evaluated by measuring .\V(z) and Vz), the input and output of the 
respective transnutting and receiving antennas. The second step necessitates factorme 
out the effects of the antennas, tn addition to elinunating the signal pollution Giaiie 
spurrous cross-coupling and clutter in the chamber. The following three measurements 


Y. 


are made with equivalent svstemn diavrams shown in [igure 
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Figure 1. The Measurement of Noise Source Impulse Response 


m Step 1. : Background measurement (no target) 


clutter 





Hi(f) = He(f) He(f) Hr(f) 
m Step 2. : Calibration target (sphere) measurement 


aa Calid. Target | 


nies ‘ 









H2(f) = Hi(f) + He(f) Hp(f) Hr(f) 
m Step 3. : Test target measurement 


target 
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= 
Hy i. 7 OY 
1 = ay 3 | 
nye! 
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H3(f) = Hi(f) + Hr(f) Hs(f) Hr(f) 
m Step 4. : Subtraction of the background effect 
Ha(f) = H2(f) - H1(f) 
Hs(f) = Hs(f) - H1(f) 
m Step 5. : Extraction of the target response (Deconvolution) 


Hs(f) = Ht (f) Hs(f) Hr(f) a Hs(f) 

Ha(f) Ht(f) Hp(f) Hr(f) Hp(f) 
For optimal compensation, use Riad's method 

Hs(f) Hs(f) Ha(f) 

Ha(f) | Hatt) (2 eee 
Therefore, 

fered Hs(f) Ha(f) 
He(f) = ~—THs(hi|2 =) COmEne yy Naa 
Hpc(f) : computed frequence response of the 


calibration target. (sphere) 


Figure 2. The Procedure of Noise Source Impulse Response Scattering Nfeasure- 


ment 


1. No target present, measure /,(2) 
2. Calibration target (sphere) present, measure /1,(2) 


Ee 1 eSt target present. ineasure /1,(7) 


Then, (2) will be subtracted from A,(2) and 4,(2) to eliminate the clutter effect. 
ha(t) = Ay(t) — fy (0) (1.1) 
fis(t) = hj (2) — Ayla) eee 


Finally, the desired impulse response is extracted from /,(¢) and /.(Q) using an opumal 
deconvolution technique known as Riad’s method [Ref. 10] which will be deseribed in 
detail in the following chapter. 

An important difficulty encountered in this thesis research was that no sampling 
device with a high enough rate for the source signal, having a bandwidth of approsi- 
fete. 13 GilzZ, is avaiable thus far. Consequently, the estimation of the crosscerre- 
lation function could not be computed using a pair of properly sampled sets of the tune 
function of cach input and output signal. Two alternate methods were developed which. 
‘for the case of ergodic random noise, permits arbitrarily slow sampling to be used. [hese 


two Nyquist independent techniques will be described in Chapter II. 


mee OVERVIEW OF TELESIS 
The objectives of this thesis were to: 
1. Develop a working software program which will acquire transient response data 
from the target and compute the estimation of the crosscorrelution function fol- 


lowed by a computation of the system inipulse respouse using a deconvolution 
technique. 


to 


Demonstrate inpulse scattering response measurements using a broad band ran- 
dom noise source, and verify the performance in comparision with another tech- 
nigue: tune-domain measurement using modified step funetion input. 

Chapter [IT will expand the theorv of noise source imipulse response measurement, 
This will exanune the basic theory of analog crosscoirelation teeliniques for acquiring 
system impulse responses and consider an appropriate derivation for the discrete version 
of the technique for the sampled signal. [nn addition, ic will quantify the estimation of 
the crosscorrelation function. 

maapter (i) describes the experimental svstem in detail. Lhis will includ: the de- 


scription of the original Digital Processing Oscilloscope (DPO) and modifications that 


were made. A description of the noise source hardware and its modified hardware will 
also be given in this chapter. In addition, the problems encountered tn this research will 
be discussed and some guidelines for the laboratory work will be summarized. 

Chapter IV contains the calibration and validation measurements for simulated 
targets. This will melude the initial quality tests of simultaneous channel sainpling, notse 
source, and crosscorrelation function that were made by sampled data using the noise 
source. The derivation of the crosscorrelation of the measurement itself 
(autocorrelation) and frequency response will be considered using PPT techniques. The 
imitial test involving a mierowave bandpass amplifier using a noise source and a step 
generator with amplifier will also be discussed. 

Chapter V will describe the electromagnetic scattering measurement. The exper- 
imental setup, scattering range and characteristics of the antenna used wil) be explained. 
Finally, the impulse response measurement of a metal sphere target will be attempted 
using the noise source. The resultant failure of this measurement will be considered. 

Chapter VI provides some conclusions about this experimental approach invelving 
noise source limpulse response measurements. Additionally, reconumendations are made 


regarding improvernents m the current system and possibilities for further research. 


il. THEORY OF NOISE SOURCE IMPULSE RESPONSE 
MEASUREMENT 


A. THEORY OF ANALOG SIGNAL CROSSCORRELATION 
Looking at the system diagram in Figure |b, the responses of transnutter antenna 
h(t) and receiver antenna /1(z) are cascaded with the parallel responses of the target /:,(7) 


and the clutter A.(7), forming an overall response, /(Z), 


A(t) = h(t) *( Add +A) * A (0) (2.1) 


The total ineasurement svstem can be represented as a srmple linear system whose re- 


sponse, }(z), due to an mnput, A(z), is given by the convolution 


F(t) = A(t) * V(2) 


as y\ 


Ee) 


= | Wa) —a)da 


a/ } 


Tt can be shown that the crosscorrelation of the input and output of a linear system 
estimates the system mnpulse response when the input haus a bandwidth that is large 
compared to the bandwidth of the system. This results from superposition as applied 
to stochastic expectation, 

Pe) = f(r) ir + 8) 
= Efh(g * X(t) A(t + 2) 
=A(1)* E(A(c).\(r + 0] 


Ky (2) = AC) * R(t) ee 


For ergodic processes, the mean values and moments can be determined by tine aver- 


Mees as well as by cisemble averaves. 


Rylt) = <\()He + 0> 


I 


| . Pte) (tc -Poneic 


=) 


Replacing the output signal }(2) by the convolution detined in Equation (2.2) results in 


Linear system 









( Fee 
AU H(t) , 
Figure 3. Crosscerrelation Measurement System Model 
RANE | M(x) | Mie Geer = adaes 
—30 0 
= | hio)| X(D)M(t + (t— 0) dedo 
0 —00 (229 
= | h(a)R{t - a) do 
0 
= A(t) * Ry(1) 
This leads to an important property of the Fourier transform relation, which is 
Sy) = AD SY) (2.6) 


known as the crosscorrelation theorem [Ref. 11]. 


B. CROSSCORRELATION MEASUREMENT 
I. Method t.: Simultaneous Dual Channel Measurement with Sampling Rate 
Greater than Nyquist Frequency. 

Three methods of performing the measurement could be used to estimate the 
crosscorrelation of the linear svstem. One possible, probably the most general, wav is 
to directly sample the input and output sequence simultaneously and store the data in 
the computer memory. This involves a delay of the input sequence, and a computation 


of the average of the product of the two sequence vectors. 


sa l reer 
Ry (72) = Ge Ss? Xe } k+n 
c=? 
K-|I—n 
Ry(n) = > AK) YK +0) (2.7) 


This scheme is illustated in Figure 4 and compared to other methods. Two conditions 
must be satisfied to use this method, 
}. Sampling must be done with a sampling rate greater than the Nyquist frequency. 


2. Sampling of two channels must be done sumultaneously. 


2. Nlethod 2.: Simultaneous Dual Channel Measurement with Arbitrary Sampling 
Rate. (Less than Nyquist Frequency) 

Method 1 could not be used because the required bandwidth of the DPO 
sampler is insufficient. An alternative method requires inanually shifang the input (or 
output) signal, being incremented by the sampling interval. Samples of the mput and 
output signal are taken with a suitable sampling rate, since cach sampled time sequence 


Meireated as the sarnple set of the ensemble space rather than a time sequence. 


A 
iy > Eta? enTk 





K~} 
eee) AKL) KEL +2) x 
a oe (2.8) 
c= 
where, 
a= AT with 2 Paes Coe ee) 
T.=LT ;' [E> (27) 


Hleie, 7 denotes the time period of delay and 7% is the sampling period. Using this 
technique, the estimation of the crosscorrelation of the delaved time point can be con- 
puted by cach measurement of input and output signal with an appropriate amount of 
delay. This method is dlustrated in Tigure 4. The remaining condition which must be 


satisfied for the crosscorrelation measurement is simultaneous sampling. 


For the case of random noise, it looks like the direct measurement method under 
the high rate sampling environment. The second method was initially chosen for this 
Tescaren 

3. Niethod 3. : Pre-subtracted Single Channel Measurement with Arbitrary 
Sampling Rate. (Less than Nyquist Frequency) 

During the laboratory work, it was shown that an unknown critical noise is 
added to the left sampling channel of the DPO which makes the second method impos- 
sible to use. This may be caused by mispositioning of the actual data m the buffer. This 
phenomena and the equivalent model are illustrated 11 Figures 5 and 6. An altermmative 
method which avoids this problem ts to subtract ¥ from WV using the built in “Add” and 
“Invert” functions in the DPO. Thus results in only one sequence of pre-subtracted data 
and is followed by computing the estimation of the mean-square value of this data se- 
quence, 

Referring to Figure 6, assume that the impulse responses of the sampler chan- 
nels of the DPO are unit impulse responses, (The effect of the sanipler resporsesioiaiae 


DPO wil be discussed in the Jater section of this chapters 
p(t) = (1) = o(?) Ret 


Then the measured sequences 1)(7) and V() can be denoted as 


A (nm) = A(t) + DLC) + NA) 2.\aap 
PF 2) a C9 a On C2) A, (2. ay 
WORE EC. 


DG) : OC bias due to the DPO vertical adjustment. 
V(t) > other noise. (thermal noise, jittem quantization 10iseiewe) 


with sunilar representations for D(z) and VQ) usig the ¥-channel. 


A. . Method 1 


N samples 
Set Kein 
x Bannan Ray(n) = » WK) YK +0) 
sete th 
Delay Samaling period 


Dual channel measurement with sampling rate greater than Nvquist frequency. 





B. Method 2 
N samples ra 
oo ! RS E > MAL) VAL +n) 
n=O 
! ’ 
Lae 
a | T,= CT 
te som = ae 
Delay Sampling period 


Duai channe! measurement with sampling rate less than Nyquist freauency. 


C. Method 3 
N samples 


1 aaa Ral) => (Ze) P-LZ)P I, 
sl 
{nt 7 = 
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single channel measurement with sampling rate less than Nyauist frequency. 


Figure 4. Comparison of The Three Estimation Techniques for The Crosscorre- 


lation 


A. Subtracted by computer B. Subtracted by DPO 
(after sampled) (before sampled) 
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Figure 5. Simultaneous Channel Sampling Test by Subtraction 
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Figure 6. <A Detailed Model of The Sampling System of the DPO 


Let a subtracted sequence be denoted as 
Z(n) = A A(AL) — YAKL + 1) (2:15) 


7 


where “n” is due to the time-shift. Then the mean-square of Z(m) is described as 


[ Zn) PP = [NAL — ¥ AL +n) 
=[{ WAL) + DAkL) + ngkL)} — {MAL +) + DKL +n) tnf(AL +a) }I 
= DC term — 22, (1) 


=[ Z(co) = Ren) 


Where we have assumed stationarity of the random process in.Y and Y. I|hereforegiim 
desired estumation of the crosscorrelation function 1s 


[[ Zloo) J? — { Zn) P| (2.14) 


A 


; | 
Rev) => 
The results of the trial measurements reveal that this estimate has a better and more 
consistent value than the result from method 2, This measurement technique haus been 


used for this research. 


C. PROPERTIES OF THE ESTIMATED CROSSCORRELAT GS 

Since the crosscorrelation function has an important role in the andiysis of linear 
systems with random inputs, an important practical problem is that of detesmiming the 
quality of the estimation of the function for experimentally observed randum processes. 
In order to evaluate the quality of this estimate it 1s necessary to deternune Cheam 
ICME a ihiee. OF Roun), since the estimated crosscorrelation 18 a radon valli 
whose preeise value depends upon the particular sample funetion being used and the 
particular set of samples taken. The mean value for Equation (2.8) and (2.1-1) can be 


computed as follows: 


l4 


k~| 
E[R,(n))= El Y MAL) MAL +1) | 


k=0 
k-l 


=> EL MAL) Y(AL + 1) ] 


k=0 


ae K R,,() 


ao 
E{Ry(1)] = Ry,() eS) 


Thus, the expected value of the estimate is the true value of the crosscorrelation (unctron 
and is an unbiased estimate of the crosscorrelation function. The vartance ts denoted 


as 


Var[ Ryn) |= E[ (Raylet) — Ry} 7 
= EY {Ry WV] -2E[ Ry ld Rela) | + EL Rg (} | 


= 5 (Ry(0)} ? — { ie yaa) }° 


—SNMAL) MAL +.) P1- Bgl) 2 


Varl R(t) |= > ELL XUL)YUL +n) } (XKLVAL +0) 3} 
Ne = = (2.16) 


Since j and & are dilferent samples in the ensemble space, it is reasonable to assume that 


they are statistically independent random variables when j #4. Hence. Equation (2.16) 


becomes 

meee GLIYGL +a} i MAL VAL +4) } p= ELL NAL) VAL +4) % Cpe 
= K.2 yt) 
= 


E{{ XUL)YUL +n) } (X(AL)Y(AL + n) } |] = ( ELYCAL) FURL +P jk 


= {R,(w} 
B 


where JA is the crosscorrelation of X? and }? while B is the square of the crosscorrelation 


for X and Y. Using this result, the variance leads to 





Var{ Ry(n) |= = KA + (KR2=K)B?] — B? 
Dele Be 
Varl R,,(n) |= = a(n) (2.17) 


Where o2,(77) ts the true variance of the product of the random variables \, andes 
An interesung phenomena takes place when the ume goes to infinity. i USaeaees 


if \(a) and F(a) are zero mean random vanables wiuien 


Rya{oo) = °F” 


Le oo) = () 


eiererore, 


2 y? 


ail Ryloo) }= = 


; 
7 (2.18) 


This means that the variance of the measured crosscorrelation approaches a constant 


value as the time separation goes to infinity. Therefore. the resultant estimated cross- 


correlation could be thought of as 


[xy ( 22) = Ry (v2) a8 : O(2) (aia) 
Where, V(2) ts the estimation error with unit vartance. [lowever, this {lect Canmiguaas 


duced by increasing the number of sample points. stnee the variance 1s normalized by 


this number. 
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D. SYSTEM IDENTIFICATION USING ESTIMATED CROSSCORRELATION 


FUNCTION. 


To measure the input and output nojse, some kind of measuring device must be 


used. Therefore it would be reasonable to think that the measured sequence results froin 


a convolution of the sequence to be measured and the svstem response of the measuring 


device. This is illustrated in Figure 7. Therefore, the measured crosscorrelation function 


1S Bayt) rather than R,(4). 


Ry (0) = | a(t +9 de 
Pere. 
i= | pla)A(t — 0) do 


= aes, 


ye) =] ae Y(e ~ 0) di 


Ae,2) 


= q(t) * ¥(r) 
Swibstituting Equations (2.21) and (2.22) into Equation (2.20) vields 


Rey (Y) = | | p(a)x(t — 9) da | g(u)y(t +e — pw) dtd 


oo 


= | ro) | | | ele) A ictal) (© <7( fl) citecte 


=| po) | | x(a)i{(4 + o — uw) + @) de g(t) duda 


p(a) | gi mee A ecitelo 


mateo Gr) 


Kye y (0) a 


DN 


a) 


wot) 


Y* 


9 fm fm 


es) 


Since the inner integral of the right expression in Equation (2.23) ts the convolution be- 


tween A,,(u) and g{u), the equation can also be denoted as follows, 


Fioure 7. 


| Noise source : 
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A Detailed System Model to be Measured 
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Rey (0) =| ne i g(t +o — W)Ry(u) dude 


—OoO 


(7 
i 
foe 
— 


7 I. [- Pla)qn(t~ w) + a) do Ry (us) du , 


lore) 
se Ryglt — H)R (Hu) dys 
Therefore, the resultant measured crosscorrelation 1s 


Rey (0) = Roglt) * Rel) 


= eC) ts R 1) a A(t) (2) 
= Ry el) * Alt) 


where Ry) iseretericulterdsailie system crasscotrclation im this thesis and has fle sane 
mole as tlhe autocorrelation tn L-quation (2.3). 

Since the measurement of the svstem crosscorrelution 1s performed without the tar- 
get, 1t should be an even function if the responses of the two sampling systems of the 
DPO are exactly the same. In fact, the responses of the two channels do not have exactlv 
the same characteristics so the test result of the measured system crosscorrelation has 
an almost even shape, as expected. This will be described in chapter 4. 


The Fourier transform of Equation (2.25) results in 


Sxy() = Sng) Sxyl(f) 
SpqV) Sx) HP) (2.26) 
= EME 


Which satisfies the crosscorrelation theorem. 
ioe set the impulse response, R, (0) must be deconvolved by R, tf). Therefore at 
ws ; *p 


least two measurements must be performed for doing time deconvolution: 


Rye xe (Mt) = Ry a(n) # R (i) (22!) 
Rey li) = Pony) = R (rt) ® fi(ae) ea 


ie power spectrums of Equations (2.27,28) are 
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Sy eK) = Spalk) SA) (2.29) 


Sy (k) = Spgl) S(t) H(1) (2.30) 


Then, the frequency response of the system 1s extracted by dividing Equation (2.30) by 
Equation (2.29). 


eee 8 
H(k) =—-———- 


I 


Se xth) 
a 2.31) 


Salk) Sk) HK) 


Spq(k) Sx(K) 


However, this frequency diviston must be compensated since the estimated cross- 
correlation function is corrupted bv noise, as described in Equation (2.19), so that the 
power spectral density is biased by the noise power. The technique used here for norse 
compensation 1s known as Riad’s Method. The “optimal estimator” prescnbed by Riad’s 
Method ike Wl ss 


where JY, represents the average spectral power density of the system crosscolrelamen 


function: 


The parameter / is denoted as the smoothing parameter and establishes the lower linut 
of the denomimiutor, hence reducing Ure notse effect on the coniputed turect 1espemee 
when the input signal spectrum is small. TPlre nnpulse response is derived by performing 


an inverse Fourier transform on the target frequency response 


h(n) = F7'{ H(k)} (2.34) 


Another rmportant point to keep m mund rs that we use an FFT algonthm which 
has quite different properties from the continuous Fourier transform. Ordinarily, when 
the FFT 1s used to estrmate the frequency response, windowing ts performed to smooth 
the result. However, applying a window interferes with the circular convolution prop- 
erty. Therefore. in this research, we will not apply windowing before performing the 
FFT algorithm. Instead of windowing, the crosscorrelation function R, (00) saliibe “Zero 
padded” to avoid the wrap-around problem which comes from the circular convolution 


property. 


Ut. BASIC EXPERIMENTAL CONFIGURATION 


A. DESCRIPTION OF DPO AND ITS MODIFICATION 

To meet the goals of this thesis work, an ultra wide band noise source and a dual 
channel sampling device with high sampling rate and precise synchronization must be 
used. A dual channel synchronous sampling Digital Processing Oscilloscope (DPO) and 
a broadband noise signal generator were used for this purpose. 

The DPO ts not an ideal sampling device since it does not sample at a rate exceeding 
the Nyquist frequency. Even though the sampling rate of the DPO is measured in MIFIz, 
it can display a periodic repetitive signal having up to 12 GHz in bandwidth by rear- 
ranging the actual low rate sampled data using a displav buffer. Hfence the signal dis- 
plaved on the oscilloscope looks like it has been sampled at over the Nyquist rate. The 
sampling technique of the Tektronix DPO sampler S-6 is tlustrated in Figure 8. The 
Tektronix S-32 pulse generator produces a fast repetiive pulse train. 7A tine wie 
output of a lumped network can be displaved at an apparent sampling rate excceding the 
Nvqurst rate by usmyg the $-52 to drive the network. This 1s shown ines 

[fowever, this property could not be extended to a non-repetitive signal Sinceuiae 
technique does not reconstruct the original shape of the actual sequence by rearranging 
the sampled data after it is taken by the sampling head. Consequently, we could not 
utilize the virtual characteristics of repetitive sampling of the DPO for this rescairchiaaaime 
broad bend noise generator has about 13 Gl1lz bandwidth which 1s far bevond the single 
shot sampling rate of the DPO. 

Using the measurement method 2. as described in Section I[.B, one of the srgnals to 
be measured ts manually shifted through the requmred sampling interval. The measure- 
nient process treats the data as samples of an ensemble space rather thai a tiem 
qucnee since our mtention rs to compute the value of the crosscorrelatron of the shilted 
time difference. Under this condition, there is no restriction of the sampling rate and the 
DO can be used. There ts, however, the need of an additronal varmable delay [iia 
on one channel to make the required tune shift. this function 1s carried out using the 


“Delav’ knob on the channel 0 sampler. 
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Figure 8. 


A. A portion of a rectangular pulse produced by the $-52 pulse generator 
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B. Modified pulse filtered by the DM ROOS bandpass amplifier 
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Figure 9. A Pulse Signal and Its Conyolved Output by a Linear Sys.om 


This delay knob has a 10 nanosecond delay range which is enough for measuring the 
crosscorrelation function for the small targets that are used. Assuming that the band- 


width of the sampled signal is less than 12 GHz, the required sampling interval 1s 


[,=7- = 40 psec (3.5) 


5 
where, according to the sampling theorem, 
fe =2B=25 Giz Bige 


To obtain this tvpe of resolution on the delay, a 10 turn potentiometer was used instead 
of the built-in one turn potentiometer on the 10 nanosecond delav knob. This 
potentiometer has a five thousand division scale on its face within the full range of its 
delay. This modification 1s illustrated in Figure 10. 

During the laboratory measurements, tt was found that the actual time delav was 
not exactly linearly proportional to the deluy reading. A potentiometer reading table 
ms tel Wscatonimmdrk the value of the reading corresponding to each time delav. 
sample potentiometer reading table is shown in Figure 11. Although the actual delay ts 
dependent upon the potentiometer reading table used with each measurement, a reading 
difference of about 3.0 represents approximately a 40 picosecond delay. 

To create a potentiometer reading table, an additional measurement using the pulse 
source was performed prior to the nolse-source measurement. After creating several such 
tables, it was revealed that the potentiometer reading has about + 3 picosecond reading 
error due to misreading and the lme-width of the displayed pulse when it 1s performed 
on 20 picosecond per division tune scule. This is illustrated in Figure 12. One factor 
Which increases this line-width is the jitter noise. This reading error increases as the tune 
scale increases and makes the potentiometer reading table next to useless. Therefore, a 
lmear increment of the potentiometer bv 1.5 1s usually acceptable if the total time win- 
dow to be measured ts over 200 picoseconds. In such a case the measurement can be 
performed using a 20 picosecond per division time scale. 

Another problem encountered during the laboratory work is a scaling dulference of 
the sampling interval of the two channcls as well as the drift of the sampling frame. A 
comparision of CY — ¥) on the 20 psec./div. and [0 nsec.. div. is shown in Figure 13. 

ooking at rigure |5, we sce that the upper part has noise with a constant cuvelope 


While the other demonstrates a noise with a consistently growing envelope. Therefore, 


it would be reasonable to think that the souree of the noise m the upper figure comes 
from the jitter noise of the two channels due to the random sampling time error while 
the other includes the jitter noise and a scaling difference of the sampling interval. One 
possible source of the scaling problem might be come from machine accuraev. So, if the 
two signals are aligned at one position using the potentiometer delay knob, so that the 
subtracted signal at the point shows a minimum envelope, then the error caused by dif- 
ferent sampling rates grows with increasurg the time diflerences froin the aligned point, 
producing larger error in the envelope. Therefore, it would be better to sample the data 
ona 20 picosecond per division (nonunal value) to mininuze problems introduced by the 
differenee of sampling rates. 

During the laboratory work it was found that the images of a pulse signal as dis- 
plaved on each of the two channels was gradually drifting from left to right with differing 
velocity. This drift problem might come [rom the thermal inconsistency of the DPO 
since it can be reduced by warning up for about one hour. However, the effect sull re- 
wnains. This problem requires completion of the experiment in as short a tine as possi- 
ee 


B. NOISE SOURCE GENERATOR HARDWARE AND IIS DEVELOR?ii ais 

A model DMR-005 broadband noise generator was used for this researel). This 
noise source includes a noise generator and two eascaded broadband ainplihers, each of 
Which has about 20 dB gain with a 13 Gllz bandwidth. The system ontput is approxt- 
mately 400 mV urrms value and the voltage histogram resembles a Gaussian distuib- 
ution. The block diagram, the frequency characteristic, and a sample recorded output 
of the noise souree are illustrated in Figures 14, 15, and 16. 

[n the laboratory it was found that a frequent change of the physical environment 
mav lead to electrical ineonsisteney resulting in output fluctuations. A little after be- 
ginning this research, one of the two amplifiers used in the notse souree failed aud it was 
sent to the manufaeturer to Ox the malfunction. After the amphiier was returned to 
NPS, an aluminum support was attached to the noise source to aveid further problems. 
The support also aets as a heat smk. Thereafter, a consistent noise output has been 


produced. 
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Figure 10. 


7) 


POTENTIOMETER READINGR ore 
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Figure Il. Potentiometer Reading Table 
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The Reading Error Due to the Line-width of Ramp Signal 


line-width = 4 psec. 
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The Effect of the Difference of Sampling Interval 


igure 13. 
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Noise Generator 


Figure I4. 
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Frequency Response of an Amplifier of the Noise Generator 
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Figure 16. A Recorded Noise Signal 
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SUMMARY OF THE PROBLEMS AND SOLUTIONS FOR THE 


CROSSCORRELATION MEASUREMENT 


The problems for the crossccrrelution measurement and the solutions described so 


far are suinmanized as follows. 


ih 


ty 


Ue 


Cay 


To produce a noise source with consistent characteristics, an aluminum heat sink 
support was attached to the noise source. 


A sampling rate lower than the Nvquist frequency and nuspositioning of the noise 
added to the channel 0 (output measurement channel) forced use of the measure- 
Niet methoues: 


The scaling difference of the sampitug period of the two channels leads to sampling 
the data on the 20 preosecond per divisiou trme scale to numinuze the eflect. 


The drift of suntpling frame requires a warm up period of at least one lrour to reach 


thermal steadv state and then to perform the nieasurement as fast as possrble. 


To perform the fastest measurement, a 1.5 reading difference of the potentiometer 
is linearly increased regardless of potentiometer reading table since the experiment 
reveals that a 1.5 reading difference of the potentioineter makes a@ 2U picoseaenm 
delay with affordable error which can be compromised with machine precision and 
visual reading error. 


‘Also, it is recomended to perform the experiment as quickly us possible te reduce 
the data transfer tune, Which 1s proportional to the number of the samples. ({t takes 
about $.5 second for trunsferrimg a set of 1024 samples of data.) 


[ollowing these guidelines, about half an hour of measurement time 1s needed for 64 


point by 2 sets of 1024 samples with an additional one hour for warnung up the svstem, 


which wul be the best condition for the crosscorrelation measurement under the ctrrent 


laboratory environment. [lowever, the drift problem stl remains since at least one-half 


hour ts needed to perform the experment. 


Finally, detailed tests revealed that the systern crosscorrelation at the peak point has 


the value of about twenty times the peak-to-peak value of the error noise. This is tus- 


trated in Figure 17. Therefore. if the value of the original crosscoriclation drops agama 


significantly, then che error noise will dominate the measured signal and will makewaie 


crosscorrelation method impossible to use. 
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error noise 


Figure 17. The Amount of the Measured Error Noise 
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IV. CALIBRATION AND VALIDATION MEASUREMENT 


A. INTIAL TEST OF SIMULTANEOUS CHANNEL SAMPLING 

The most important and critical role of the DPO in this research 1s in simultaneously 
sampling a pair of time functions which are probed at the input and the output port of 
a linear system. A useful validation test for the simultaneous sampling characteristic is 
to subtract the signal from one channel (VY: input secjuence in this thesis) from the other 
channel (2: output sequence). This will be displaved on the DPO screen by using the 
burt in “Add” and “Inverse” functions of the DPQ. Using the manual delay knob. the 
image of the subtracted signal can be mininuzed using the simultaneous sampling posi- 
tion. The potentiometer reading at this point (1t must be the 7 =U point of the system 
crosscorrelation, R,. (2) has a value of about 200 indicating that the sampling channel 
Q (output channel) has a tune delay of about 2.7 nanoseconds less than that of the 
channel |. This is illustrated in Figure [8 with a comparison with the original tnjamiies 


qpuenece, 


Looking at the Iigure [8, the left part shows a noise with a consistent envelope of 


small variance. The source of the noise on the left figure comes from the jitter noise of 


the two channels produced by random sampling time errors. Assuming that the enyel- 
ope of the jitter noise is proportional to the time derivative of the source signal itself, the 
noise source signal must have more variance of tts jitter noise vis-a-vis the ranip source 
signal. This results since the time derivative of the ramp signal produced by DPO is 
much smaller than that of the noise source indicating that the noise signal is more vul- 


nerable to the jitter noise than the ramp signal. This property is verified in Figure 19. 
BG. CROSSCORRELATION OF THE DPO Sys 

‘As deseribed in Section I1.D, the crosscorrelation between input and output 1s con- 
volved by the crosscorrelation between the impulse response of the two sumpling heads. 
Tine ciscrels version of cas anals 
P = a 
Key (1) = Kaglit) * Kyl) 


= R(a) * Ria) * h(n) (4.1) 


= Re (1) = ft( 42) 


A. Channel | - Channel 0 ere banner 





SCALE: @.1 VOLT SCALE: @-E VOLT 
MEAN 2: G.a3S8273625125 MEAN : BS. G1G1g?re79eg3 
UAR : 7.428921318E-4 UAR + B.0106785155699 


Figure 18. Simultaneous Channel Sampling Characteristics by Subtraction 
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Jitter Noise of the Noise and the Ramp Signal by Subtraction 
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igure 19, 
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where, 


AN 


R, «() = R,,(n) * R(n) (4.2) 


p 


mere R, (2) is the estimation of the system crosscorrelation and ts convolved with the 
crosscorrelation of the sampling head and the autocorrelation of source noise itself. To 
approximately measure the system crosscorrelation, we simply feed the same input signal 
to each sampling head, as illustrated in Figure 20. 

In the ideal case. each sampler svstem impulse response, p(1) and q(t), must have the 
same shape, and at best the unit impulse response. Although these two inipulse re- 
sponses are not exaetly identical, they will mn reality be very close to each other. Conse- 
quently, the shape of the system erosscorrelation R, (7) IS (ClOSe Ho DUl ier CaCl dtl 
even function, while the autocorrelation of the source noise R(n) must be even. 

The initial ineasurements were performed using a 2.5 picosecond sampling period to 
get a more detailed time shape of the system crosscorrelation. Subsequent measurements 
Were performed using a 10 picosecond sampling interval to see the entire shape of the 
system crosscorrelation. Figure 2t and 22 show the test results. As expected, the shape 
of the crosscorrelation is almost even, which means the two sainpling channels have 


Bar characteristics. 


C. DECONVOLUTIGN VALIDATION TEST 
I. Test with Noise Source 

An Avantek model SAS3-2954 solid-state inicrowave amplifier was used during 
this validation test. This amplifier has an average gain of 42.5 dB over the band from 2 
to 6 GHz and provides a 3 dB bandwidth of about § Gliz. The spectral characteristic 
of the Avantek amplifier is shown in Figure 23. This amplifier is linked to two 20 dB 
attenuators, since it produces moderate power, to prevent anv possible damage to the 
sampling heads which have an operating range of 5 volts peak. Figure 244 illustrates the 
basic setup for this test. 

To extract the mnpulse response of the Avantck amplifier (2-port), the following 
three steps were followed: 

1, Test without the 2-port, measure the system crosscorrelation. 


2. Test with the 2-port, measure the crosscorrelation of the 2-port. 


3. Perform deconvolution computation. 


oy) 


The system crosscorrelation measurement was performed using 64 time points 
having a 40 picosecond increment. This gives about a 2.5 nanosecond window, which ts 
almost three times longer than the effective period of the system crosscorrelation. The 
amount of the delay time interval (40 picoseconds) ts sufficrently small in this case smece 
this gives a Nyquist frequency exceeding the bandpass frequencies of the Avantek am- 
phifier. A plot of the measured system crosscorrelation ts shown in Figure 25. 

The bandpass amplifier crosscorrelation measurement was performed using the 
same tine points and delay interval. Further, an initial delay was given to the channel 
Qin the amount of 400 picoseconds to compensate for the effective delav through the 
Avantek amplifier. The resultant measured crosscorrelation is shown in Figure 20. 

Finally, the computation of deconvolution was performed using Riad’s optimal 
compensation technique. As deseribed in Section [1.D, this deconvolution method 
serves to limit the lower bound of the denominator and thus reduces the noise effect. 
Several coniputations were performed with different smoothing factors and finally the 
cuse Of 0.1 was selected for the best solution. After the computation of deconvolution, 
the computed impulse response was moved 32 points to the right to show the best 
graphic interpretation. The result is illustrated in Figure 27, 

2. Test with Modified Pulse Source 

The same three measurement steps were performed usiig the tine dona 
method described by McDaniel [Ref. 9]. The output ramp signal from the S-32 pulse 
generator was filtered by one of the amplifiers used with the norse gencrator to produce 
a roughly similar power spectral density at the input. This source signal has enough 
Innate Zero pacding to elinunate the wrap-around problem. The test was fairly straight 
forward and the results are shown in Figures 28, 29, and 30. 

The comparison of the two methods is given m Figure 31. These twosresmi 
have almost the same shape, as is expected, aside from the time scaling difference. Ihe 
crosscorrelation method uses about a 4O picosecond time interval while the time domain 
inecthod samples with a period of 39 picosecond, which 1s fixed by the DPO spaces. 
Another fact is that the result from the crosscorrelation method has more error gee 
power than that of the pulse method. Therefore, the performance of the crosscorrelation 
method using the current DPO appeurs to produce lower quality results than that of the 


pulse method. 
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Figure 20. DPO Svstem Crosscorrelation Measurement Setup 
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Figure 21. DPO System Crosscorrelation Function 
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Figure 22. DPO System Crosscorrelation Function and Power Spectral Density 
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Frequency Response of the Avantek Amplifier 
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Figure 24. Noise Deconvolution Validation System Setup 
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Figure 25. The System Crosscorrelation, Crosscorrelation Method 
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Figure 26. The Crosscorrelation of Avantek Amplifier, Crosscorrelation Nlethod 


47 


HUMBER OF SHIiPbS POL 


6&4 


SPEC 


FEST WITH 
TRRGET 


Si Glgeeec) S 
Ri on aan 
tes ee 
MAA. 3 UY 
fg AR) 
BYitec 3 1 


SMP. < Souter 
cc KES. : 6. 23 biee 
Eee WIH. 3 12 34G shies 
MAK. 2 4.105 
SHUUTH EG eee 
ROTHITE FS Sees 


Figure 27. The System Response of Avantek Amplifier, Crosscorrelation Nlethod 
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Figure 28. The Modified Pulse Input Signal, Time Domain Method 
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Fieure 29. The Output Signal of the Avantek Amplifier, Time Domain Method 
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Figure 30. The System Response of Avantek Amplifier, Time Domain Method 
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Figure 31. The Comparison of the Two Methods 


V. ELECTROMAGNETIC SCATTERING MEASUREMENT 


A. DESCRIPTION OF THE SCATTERING RANGE 
To derive the characteristics of the noise source impulse scattering response meas- 
urement, a physical environment and equivalent model for the scattering range should 
be developed. Since the scattering range was already described in detail by Mariatequi 
feet. > |, and NicDaniel [Ref. 9], only a bricf description will be given in this sectton. 
The general physical environment of the scattering range ts illustrated in i igure ta. 
As shown in the dlustration, the svstem can be divided into 5 parts. 
1. Anechoic chainber 
2. Transnutting and receiving antennas 
3. Broadband analog noise generator 


4. Dual channel sampung device (DPO) 


tay 


Signal processor and controller (Tektronix 4052. desktop computer) 


— 


The purpose of the anechoic chamber is to shicld the incident and scattered 
waveforms from external noise and interference effects as well as multipath contam- 
mation. The physical dimensions of the chamber are 

eeLoncitudinal lengtl: : 20 ft. 
ee Laicral width: 10 [t. 


merit : 10 ft. 


Ca 


To provide shielding from atmospheric and man-made noise, the chamber 1s 
internally covered with alunnnum panels which: are earth grounded. A special absorber 
material manufactured by Pantec, a division of Emerson Electric, is attached to the 
alununum shield to absorb the electromagnetic radiation. This absorber is ridged along 
Mmemenett: of the chamber to cure the incident wave to the back wall where the most 
of the incident and guided clectremagnetic wave is absorbed by 18 inch pyramidal cones. 
Eteht-inch pyramidal cones are attached to the front wall, where the antennas are 
mounted. These absorb back-scattered radiation and prevent reradiation of the 
electromagnetic wave. This material is designed to provide back-scattering attenuation 
ofa 500 MEIz signal by about 30 dB below that of a flat metal plane. Reflection coctfi- 


cients are increased to about -12 dB at 1 Gllz for near grazing incidence on the side- 
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walls, floor and ceiling. This leads to a realistic low frequency limtation of about 1 GHz 
due to the resultant multi-path interference. 

A styrofoam column is used as the support for the various targets in the chamber. 
This material has a relative permittivity of about 1.1 at 3 GIIz so that it gives neghgible 
effect (reflection, refraction) to the transmitted wave but provides a stable support for 
the targets. 

Two double-ridged horn type antennas are used for the transmitting and receiving 
antennas. They have a usable bandwidth of | to 12.4 GHz and have a relatively flat gain 
over the band. 

The broadband analog noise generator, DPO, and the Tektromx 4052. desk top 


computer were alreadwadescribed in Chapter 111: 


GB. DERIVATION OF NOISE SOURCE IMPULSE SCATTERRING RESP ORE 
MEASUREMENT FORMULA 

Since the goal of this research is to validate the crasscorrelation measurement DV 
comparing the impulse scattering response result with that from the trme domain meas- 
LITeMNENL, a Sinimle sphere Wwasrenocen [or the tarccer 

An equivalent system diagrain for the scattering range is shown in Figure G2 seein 
the rmpulse responses are 


lL. At) : Target 


2. A(t) + Clutter and antenna cross-couphng 
oa Att) = Transinit-antenna 
4. Ate)  : Receiving antenna 


and the stgnals are 


tr) : Transmutted (source) signal produced by DPO 
Bey Oe : Received (back scattered) signal 


3. (kL) + Sampled transmutted sequenee 


4. u(AF +): Sampled received sequence 


(an 


z{n) 7 Stored secuemee (lo 1) 
6, Ry (7) : Crosscorrelation output 
Compared with the discussion in Section [V.C, the system impulse response in the 


anechoic chamber is more complicated than a simple 2-port network responsenmaane 


impulse response of the anechore chamber ts 
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Figure 32. A System Diagram of the Scattering Range 


A(t) = Ado) * (s(t) + ACD) * A(t) 
= (A(t) * A(t) * A(t) + Ade) (5.1) 
= A(t) * ast) + hy) 


This equation suggests the need of at least three sets of measurements for extracting the 
target response /,(r). Using our method, the three crosscorrelations that must be meas- 
ured are 
Rey f Pe lene tee (a 
=| Reg) © Re ee) (See) 
= (p(t) ‘ Aft) 


Ry y (io = (0) * [Ag(0) + hel) | (5.3) 
Rey (i)3 =P(1)*(A()+ h(n) | (3.4) 
Mere. 

OED (ar aa isl) ae Hey ar 9 (3:59 


liowever, experiments revealed the nced of one more background measurement for bet- 
ter subtraction of the background noise. This is because of the use of a relatively short 
time Window compared with the difference of the propagation delays of different targets, 
(the calibration target, and real target). Therefore, instead of Equation (5.2) to (3-4), 


we use Equation (5.6) to (3.9). 


al = tnt) * h(n), (5.6) 
Ry (13 = b(n) * ioe (Sa) 
Rey (it)y = ht) #[hyln) + hang} (3.8) 
Re. ( 2) = b(n) *[ A(n) + hay | (5.9) 
where 


P(t) = Ryglit) * Rx(n) * h(n) * h,{n) (5.10) 
h(n); cS h(n — m), (G98) 


Here, #7 is the time delay due to the different measurement tine origin caused by the 
physical shape of each target. So, the best delav could be determined by the dillerence 
of the transit time of the leading edge of the scattered electromagnetic wave for each 


target. 


m= Integer { x } 
~ Ole 

In this equation, D denotes the radial difference of the two spheres and c represents the 
free space velocitv of an electromagnetic wave. This is illustrated in Figure 33 and an 
Smiore Of this elfect is shown tn Pigure 34 using an | sq. ft. copper sheet target (or 
Petter rellectioin. 

lite ext thing to do alter the four measurements of the crosscorrelation is the ex- 
traction of the target response /,(7). This is performed in tvo steps. 

ieec iutter subtraction. 


2. Crosscorrelation deconvolution using Riad’s method. 


Pat first, subtraction of the clutter effect was performed. Let the subtracted result be 


denoted as 


R31 (77) = Ryy (a = Rey {12 


(5.12 
= (12) * hn) 
R34 (a) a Ryy (13 a Rey (3 
fo: 13) 


The following is the deconvolution for extracting the response of the target. The 
Rey pot of the deconvolution is the use of the calibration target which has a known 
computed transfer function. The Mie-series program prepared by Professor Morgan 


computes the transfer functron of the sphere [Rel. 9]. 





Figure 33. 
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1624 


Let /7,.(f) represent the computed transfer function of the 8 inch sphere. For an ideal 


case, the transfer function of the target could be computed as 


Sai) 


AW = SH 


H,<(f) (5.14) 
Where S,,(//) and S,,(/) are the Fourier transforms of the crosscorrelations R,,(7) and 
R,,(). Since we are dealing with the estimated crosscorrelations, Riad’s method, as de- 


scribed in Section [.D, ts applied. 


’ Sue: 
H.{k) = eee 11,k) (5.15) 


[Sei 


where, 


Qn 


Finally, the smoothed estimation of the mipulse response will be obtained from the tn- 


verse Fourier transform of Equation (5.15). 


C. ELECTROMAGNETIC SCATTERING GROSSCORKE Rite 
NWMEASUREN LENG 

First, a 12 inch calibration target measurement of R,, (1); was performed. To esti- 
mate a suitable time windew, a pulse signal was first used as the excitation in the 
anechorc chamber. The measured backscattered signal contains several constituents: di- 
rect antenna couplmg; chamber clutter and target scattering which 1s modired by the 
antenna responses and multipath. The target pedestal is located approximately $ [cet 
from the horn antennas. 

An additional 7 feet of cable delay line was linked between tlre transinitter power 
splitter and the channel 0 sumpler to show the two signals m the same time window. At 
first, a directional coupler was selected instead of the power splitter stnce the directional 
coupler transnits most of the input sigual. Llowever, it Was found that the directronal 
coupler provided a significant attenuation in its high frequeney band. The result of the 


See’ 


alignment test 1s illustrated in Figure 35. The leading edges of the two signals are aligned 
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when the potentiometer reading indicates 195. The time width of the sphere response 
Was about 2 nanoseconds. 

Next, the noise source was applied to measure the electromagnatic scattering cross- 
correlation measurement. Since the pulse response 1s limited to about 7 Gliz [Pig. 23.], 
a set of 64 samples with a sampling mterval of 40 picoseconds was chosen to cover a 2 
nanosecond time interval. This produces a Nyquist frequency of 12.7 GFElz. 

Unfortunately, the measured value of the crosscorrelation was significantly cor- 
rupted by the noise and system errors so that it was impossible to use. Since then, the 
Bien was focused on finding the source of the errors and on the wavs to avoid it. 
Meorvever, tlic answer is that, at this tyme, we could not reduce the error sulliciently to 
makes the estimation useful. In fact, the fight against the error noise for the signal noise 
source spans the whole history of this research. The situation tn the seatlering meas- 
urement offers signifcantly less output signal strength than was observed in the initial 
tests using lumped filters. This degradation of signal appears to have exceeded the limits 
of fidelity needed to demonstrate any viability of the measurement. Jn short. the results 
appeared as “garbage” that could not be reproduced on anv two measurenients. The 
Sichial was essentially buried in the noise. 

In addition to the description about the noise in Chapter I], the sampled sequence 


of the two channels are modeled as 


M(t) = N(0) * p(n) + jlrr) + Oeln) + dls) + Nel) + Dyn) (5.17) 
ee) = (1) * g(a) +1) + Od) + aya) + 5,(n) + u,(n) + Ni(a) + Dyas) feelc) 
ere, 


1. A(z): Input signal 


tr 


p(n): Impulse response of the input channel sampler 
3. Ja): Jitter noise of the input channel 

4. O07): Quanuization noise of the input channel 

5. d(n): Drift (nonstationary) noise of the input channel 
fee 7): Otlier noise (Thermal ioise etc.) 

WD {nj}: DC bias due to the DPO vertical alignment 

pee (7): Output sicitdl 

9. g(t): Lmpulse response of the output channel sampler 


10. y,(77): Jitter notse of the output channel 
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11. O72): Quantization noise of the output channel 

12. d(m): Drift (nontationarv) noise of the output channel 

13. s,(2): Sampling interval scaling noise of the output channel 
4. w(m): Mispositioning noise 

15, N(72): Other noise (Thermal noise etc.) 


16. Dv): DC bias due to the DPO vertical alignment 


[t may be surmised that the poor performance 1s probably caused from the nonstation- 
ary characteristics of the DPO, as sunwmnarized in Section IH.C. This inelides teams 
noise and quantization noise. The estimation error for a stationary process should be 
reduced bv increasing the number of sample points. As shown in Figure 36, the error 
noise 1S not linearly reduced by increasing the number of the samplcs. This reveals the 
existence of other sources of the error such as nonstationary noise (caused by the drift 
problem) which may give a sudden level of noise that can not be reduced by fiieiaia eae 
UO 

Figure 36 reveals that the clectromagnectic scattering crosseorrelation measurement 
for the 12 inch sphere is significantly corrupted by the error noise although the result 
came {from the computation of 8192 samples. Therefore, 1t might be almost impossible 
to perform deconvolution under such a low SNK environment unless other ways can be 
found to increase the SNR of the estimated value. One such way. which may soieumae 
problem, 1s the use of a new generation of sampling devices, such as the Hewlett Packerd 
Nfodel 54120 digitizing oscilloscope. This DPO has a sampling rate up to 20 GJ1z and 
offers computer controlled time-offsets of its 4 channels. A follow-on investigation will 


consider the application of the HP-DIO to this research. 
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Figure 36. The Fluctuation of the Estimated Crosscorrelation Due to the Error 


Noise 
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A. SUNIMARY 
The main effort described in this thesis 1s to demonstrate the viubilitv of unpulse 
response measurements using a broadband random noise source. An application of this 
system is radar target detection and tracking, where it would offer the following potential 
advantages: 
feeleOVseost Noise sources are available. 


2. Deception in radar applications, (jaminer or natural noise appearance). 


Mie limtstep cf the pesearcn Was to Gevelon the theory of noise source tmpulse re- 
sponse measurements. An analog version of the crosscorrclation was derived followed 
by development of the three crosscorrelation measurement methods. The third incas- 
urement method (pre-subtracted single channel measurement with arbitrary sampling 
Sie) was selected alter much eflort m trying to use method 2. The bias and vurianec test 
Situc estimated crosscorrelation lollowed. [t was shown that the estimation 1s unbiased, 
and consistent. finally, tle effect of the sampler on the signal as well as the application 
of an optimal deconvolution technique were considered. 

Sinceientesec was devotcd to the experimental configuration. The churacteristics 
Seecie iektronix DPO and noise generator were investigated. [t seems that the 
Tektronix DPO is not the ideal device for this research since it exhibits many problems 
for these measurements. Some of these are the drift problem, the relatively large amount 
of jitter noise (about 1/5 of the signal strength), a sampling interval scaling problem, ete. 
An cffort to minimize this error noise was attempted. Detailed tests were performed. and 
fects revedicd that the system crosscorrelation included errors in its peak-to-peak value 
ef about 1 20 of its maximum signal power. 

The calibration and validation test were then performed. An .Avantek solid-state 
nucrowave amplifter was selected as the test item. The crosscoirelation measurements 
of the sarmphng system and the Avantck amplifier were performed. A smocthed fie- 
quency response of the Avantek amplhiicr was computed using the Riad’s method. 
Mien, the impulse response of the Avantck amplifier was obtained using the inverse 
lourier transform. ‘This was compared to results obtained via a direct pulse excitation 


measurcinent. 
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Finallv, an electromagnetic scattering measurement was attemped after considering 
the requisite mathematical foundation. Four measurements were needed in practice al- 
though the theory requires three measurements. The scattering measurement provided 
insufficient fidelitv due to hardware linutations of the Tektronix DPO. Major sources 


of error in the current system include nonstationary drift, jrtter,and mispositioning. 


B. FUTURE CONSIDERATIONS 

The main difficulties encountered in this research resulted from the use ofa less than 
ideal sampling device that could not sample sulficrently fast to satisfv the Nvquist crite- 
rion. Turther, it has been revenled that this sampling device has many problems as de- 
scribed in the previous section. 

The significant payoff of this research mandates that it be continued with vigorous 
attempts to alleviate the hardware deficiencies eicountered in this initial effort. The use 
of a new generation digital processing oscilloscope by Hewlett Packard is a step m the 
right direction since this device allows computer control of the temporal offset of each 
channel with an accurate resolution of 10 picoseconds. The essentially drift-free pei 
formance, coupled with automated time-olfsets mn the crosscorrelation data acquisition, 
wil allow large ensemble averages to be employed with much shorter measurement tines 


than have been the case here. 
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